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ABSTRACT 

Recent observations of the environments of gamma-ray bursts (GRBs) favour massive 
stars as their progenitors, which are likely to be surrounded by gas and dust. The 
visibility of the optical and UV emission of a GRB are expected to depend on the 
characteristics of both the dust and the GRB emission itself. A reasonable distribution 
of surrounding dust is capable of absorbing all the optical and UV emission of the 
optical flash and afterglow of a GRB, unless the optical flash has a peak isotropic 
luminosity L pca k > 10 49 ergs _1 . This means that dark bursts should exist and these 
bursts will have to be studied at infrared rather than optical wavelengths. In this 
paper details will be given about the infrared GRB dust emission. The reprocessed 
dust emission peaks at a rest-frame wavelength of about 8 /xm. Forthcoming space 
telescopes, in particular the IRAC camera aboard the Space Infrared Telescope Facility, 
could detect this emission out to a redshift of about 2. However, an accurate position 
of the GRB afterglow must be provided for this emission to be identified, because the 
light curve of the reprocessed dust emission does not vary on time-scales less than 
several years. 
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1 INTRODUCTION 

The first gamma-ray bursts (GRBs), short and intense 
bursts of 0.1 — 1 MeV photons, were discovered at the end of 
the 1960s (Klebesadel, Strong & Olsen 1973). In the 1990s 
the BATSE detector on board the Compton Gamma-Ray 
Observatory^ observed over 2500 GRBs and showed that 
GRBs were distributed isotropically over the sky. A sus- 
pected conclusion was that GRBs were of extragalactic ori- 
gin at cosmological distances (see for example Paczynski 
1995). This conclusion was confirmed in 1997 using the Bep- 
poSAX satellite by the discovery of a well-located transient 
X-ray counterpart (Costa et al. 1997), and several hours 
later its association with a fading optical point source (Groot 
et al. 1997) at a redshift of 0.695 (Sahu et al. 1997). Many 
more afterglows have now been discovered and more red- 
shifts have been measured (Greiner 2001). Based on the ob- 
served gamma-ray flux densities, this implies that the total 



see ittp:/ /cossc. gsfc.nasa.gov/cossc/cgro/batse_src. html 



energy radiated by a GRB is of order 10 to 10 erg, if the 
emission is isotropic, and is released in a few seconds. 

The existence of multi-wavelength afterglows was pre- 
dicted by models that describe GRBs as opaque radiation- 
ionized fireballs, generating a relativistic blast wave (e.g. 
Meszaros & Rees 1993, 1997; Katz 1994); see Piran (1999) 
for an extensive review of the fireball model. The most 
promising candidate objects that generate GRBs are bi- 
nary neutron star mergers (Piran 1999) and hypernovae or 
collapsars (MacFadyen & Woosley 1999; MacFadyen 2000 
and references therein). Observations of the environment in 
which GRBs occur circumstantially favour massive stars as 
progenitors of GRBs. For example, multiband photometric 
imaging of four GRB host galaxies show their that spec- 
tral energy distributions (SEDs) are best fitted by starburst 
spectra (Sokolov et al. 2001). Also, there are some indica- 
tions from Hubble Space Telescope observations that GRB 
afterglows are associated with bright blue star-forming re- 
gions in their host galaxies (for example Bloom et al. 1999b; 
Fruchter et al. 1999; Bloom, Kulkarni & Djorgovski 2001). 
The Infrared Space Observatory observed the location of a 
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single GRB, GRB970508, both weeks and months after the 
burst and detected non-transient emission at 60 /im, which 
was interpreted as an indication of the presence of an ultra- 
luminous infrared galaxy (ULIRG) (Hanlon et al. 2000). If 
this ULIRG is the host of the GRB, then the estimated star 
formation rate must be around 200MQyr -1 , likely placing 
this GRB in a very powerful star forming region. Another 
potentially important discovery, which favours massive stars 
as progenitors of GRBs, is evidence of a possible supernova 
component in the afterglow of GRB970228 (Reichart 1999) 
and GRB980326 (Bloom et al. 1999a). However, an alter- 
native explanation for this, in which scattering of a prompt 
optical flash by diffuse dust surrounding the GRB explains 
the observations, has been proposed by Esin & Blandford 
(2000). 

If the progenitors of GRBs are indeed located in star- 
forming regions, then the GRB is likely to be surrounded by 
gas and dust. This can influence the spectrum of a GRB, 
and there have been several recent studies of the effects. 
Meszaros & Gruzinov (2000) showed that a dense environ- 
ment should modify the X-ray afterglow of a GRB, flat- 
tening the X-ray light curve from hours to days after the 
burst, before steepening the decay at later times. Waxman 
& Draine (2000) investigated dust sublimation by GRBs 
and its implications. They found that a prompt optical flash 
could destroy up to 10 7 Mq of dust out to ~ 10 pc around a 
GRB. They also calculated that the near-infrared luminos- 
ity from dust heated by the optical flash could be as great as 
10 41 ergs" 1 . Perna & A guirre (2000) proposed to use GRB 
afterglows as probes of galactic and intergalactic dust. Re- 
cently, Galama & Wijers (2001) found high X-ray absorption 
column densities, but low optical extinctions for their sam- 
ple of GRB afterglows. They inferred that most GRBs are 
embedded in large molecular clouds, and are therefore likely 
to be produced by the death of short-lived massive stars, 
but that the surrounding dust is destroyed by hard radia- 
tion from the GRB, leaving the afterglow unextinguished. 

In this article, the influence of surrounding dust on the 
visibility of the afterglow and optical flash of a GRB is in- 
vestigated in order to study if such luminous events could be 
detected at high redshift, even if they occur in dense dusty 
environments. The interaction of homogeneously distributed 
dust around a GRB with the optical and UV emission of the 
burst afterglow is modelled. The unextinguished light curve 
of a GRB is described in section |i| A dust model consisting 
of carbon grains with a range of sizes is used, which takes 
account of the evolution of opacity, dust sublimation and 
temperature; see section ^| The model is used to estimate 
both the visibility of the optical/UV afterglow for different 
densities and profiles of the surrounding dust (section |^), 
and the properties of the thermal emission of this dust (sec- 
tion ||). Finally the flux densities of the dust emission are 
compared with the sensitivity of future infrared detectors 
(section |^) . 



2 GRB OPTICAL LIGHT CURVES 

The optical light curve of a typical GRB consists of two 
parts: a prompt optical flash coincident with the gamma- 
ray emission and a late-time afterglow. The afterglow can 
be explained by synchrotron emission from a decelerating 
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Figure 1. B— band flux densities 1 Mpc away from a GRB 
with a light curve described by the model of Sari et al. (1998) 
used here. An afterglow with a Gaussian optical flash with 
^peak = 10 49 ergs -1 in the energy range hv = 1 — 7.5eV is 
shown by the solid curve; one with a Gaussian optical flash with 
^pcak = 10 48 ergs -1 is shown by the dashed curve; an afterglow 
without an optical flash is shown by the dot-dashed curve; and an 
afterglow with a top-hat optical flash with -L pca k = 10 48 ergs -1 
is shown by the dotted curve. 

relativistic shell in collision with an external medium: see 
Sari, Piran & Narayan (1998) for an analytical description 
of the light curve or Piran (1999) for a more detailed deriva- 
tion. 

We assume that the peak frequency of the optical flash 
falls in the optical waveband (Sari & Piran 1999), and that 
the time profile of the optical flash is represented as either 
a Gaussian or a top-hat function. Both the peak luminosity 
and the duration of the flash are left as free parameters, but 
typical values of the peak luminosity are L v = 0.5 — 5 x 
10 33 ergs -1 Hz -1 at a frequency v = 1.4 X 10 15 Hz, with a 
SED L v oc v" 1 / 2 and a duration of order 10 s (Akerlof et al. 
1999; Waxman & Draine 2000). 

Waxman & Draine (2000) showed that photons in the 
energy range 1 — 7.5 eV will dominate the heating of the 
dust, as in dense regions (typically with hydrogen column 
densities nu > 10 2 cm -3 ) photons with energy hv > 13.6 eV 
will photoionize hydrogen very efficiently, while 7.5—13.6 eV 
photons can be efficiently absorbed by H2, by either vibra- 
tional excitation or photodissociation (see also Draine 2000). 
The number of photons at hv > 7.5 eV that are available 
for dust heating will thus probably be negligible. The lower 
limit of 1 eV is taken because both the fraction of the energy 
radiated by the GRB below this energy is small and the ab- 
sorption efficiency of dust decreases rapidly at wavelengths 
longer than about 1 /j,m (= 1.24 eV) (Draine & Lee 1984). 

The following parameters are used to describe the after- 
glow, based on values derived from observations (including 
Galama et al. 1999a; Lamb, Castander & Reichart 1999; 
Huang, Dai & Lu 2000; Greiner 2001). The energy of the 
shell E = 5 x 10 52 erg, the fraction of the shell energy that 
goes into the kinetic energy of electrons e c = 0.15, the ratio 
of the magnetic field energy density to the total thermal en- 
ergy £b = 0.1, the surrounding density ni = 5xl0 3 cm -3 and 
the coefficient of the power-law distribution of electron ener- 
gies in the GRB shell p = 2.5. The hydrodynamic evolution 
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of the shock is assumed to be adiabatic, because radiative 
evolution only applies if e c ~ 1. 

In Fig. hi four different optical light curves of GRBs, 
used in our model, are plotted. 



Now, using equation (gj) and including the sublimation, 
equation (^]), the out-of-equilibrium temperature Td of a 
dust grain can be calculated using 



3 RADIATION ABSORPTION AND EMISSION 
BY DUST GRAINS 

3.1 Dust temperature 

The temperature of dust determines its emission spectrum. 
Because the GRB afterglow is very luminous, the influence of 
heating by both reprocessed dust emission from within the 
dust cloud and background radiation from the interstellar 
radiation field of the galaxy can generally be neglected. The 
dust grains are assumed to be spherical with a radius a. A 
source of radiation with a luminosity L = f L v dv is at a 
distance r from a dust grain. If the temperature of a dust 
grain is less than the temperature at which sublimation is 
a more efficient energy loss mechanism than radiation, then 
the dust temperature T<j can be calculated using 



i L 
4nr 2 



3 (a) = 4ira Q e 



(1) 



(all the equations in this section are adapted from Drainc 
& Lee 1984, Siebenmorgen, Kriigel & Mathis 1992, Evans 
1993 and Waxman & Draine 2000). Q a bs is the weighted 
effective absorption efficiency, a is the Stefan-Boltzmann 
constant, Cv is the specific heat capacity and Qom is the 
Planck-averaged emission efficiency, 
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in which r is the optical depth between the source and the 
dust grain (described in more detail in section 3.3). Equa- 



tion (H) has the same structure as equation (8) of Waxman & 
Draine (2000), but with a term included for the heat capac- 
ity, which can be important at temperatures below the sub- 
limation temperature, where thermal emission is the most 
important cooling mechanism. 

Equations (^) and ^ are coupled differential equations 
for the evolution of the temperature Td(t) and radius a(t) of 
a dust grain. 

The period of time during which a grain is entirely sub- 
limated after the arrival of a GRB afterglow is defined as 
the destruction time and is, in general, a function of both 
grain size and distance to the source: ides = td cs (a,r). 



3.2 Dust luminosity 

The luminosity of a single dust grain is 



x(a,T) = 



J B v {T)Q v ^ hs {a)dv 



fB„{T)dv 

The weighted effective absorption efficiency 
Jo L v Q Vi3 b s dv 



Qa 



(2) 



(3) 



depends on the spectrum of the source. 

When subjected to very intense radiation from the opti- 
cal flash of a GRB, the temperature of dust can rise to levels 
where cooling of the grains by sublimation is the dominant 
cooling process. 

The cooling rate due to sublimation of a spherical dust 
grain per unit time, -E su b, 



2 p „cla 
47ra — B — 
m dt 



(4) 



can be described in terms of the reduction in the grain ra- 
dius, da/dt, where p is the mass density of the grain, m the 
mean atomic mass and B is the chemical binding energy per 
atom within a grain. An approximation to the sublimation 
rate, da/dt, is given by 



da 
dt 



Poe 
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Po and To are constants, which depend on the properties 
of the grain material. 



L d {a, t) = 4na{t) 2 Q cln [a(t) , T(t)] aT{t) 4 



(7) 



assuming the grain radiates as a modified blackbody with 
an emission efficiency Q cul (equation ^) . 

The flux density received a distance r from a grain, 
that is radiated in a frequency band between ^ m i n and f max , 
can now be calculated. Defining x mln — hv m i n /kT(t) and 
i m „ = hv max /kT(t), the flux at a distance r from a single 
grain is approximately 



F{a,t) 



1.9xl0~ 2 L A (a,t) 
4nr 2 



■ dx , 



(8) 



assuming a reasonable emissivity index j3 of 1.5, where 
e v oc v 13 . To calculate the evolution of the flux of a shell 
with radius R consisting of grains with size a, the flux of 
all the grains must be added, taking into account the time 
delay due to the differences in distance from the individual 
grains in the shell to the observer. A photon emitted by a 
grain at an angle 6 from the line of sight has to travel an 
extra distance of R(l — cos 8) as compared with a photon 
emitted from the same shell at the same time, but directly 
along the line of sight in order to be detected simultane- 
ously by the same observer. If we assume that the distance 
to the observer is much greater than the radius of the shell 
(r ^> R), and the number density of grains is n, then the 
total flux received by the observer is the integral over all 
shells: 
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F toX {t)= f [ F(R,a,t)dRda = 

J a J R 

[[[ 2nR 2 n(R,a)sm9F[a,t+-(cos9-l)]d6dRda. 

J a J R J 6 C 

(9) 

Equation (^) can only be used for calculating the dust 
emission in the case of a spherical distribution of dust around 
a GRB, the emission from which is itself isotropic. If, for ex- 
ample, the afterglow emission of the GRB is beamed, then 
(roughly) only the dust inside the opening angle 9q of the 
beam will be heated. In that case, the emission will be re- 
duced by ~ #o/4, and the time profile of the dust radiation 
will change, depending on both the opening angle and in- 
clination of the beam to the line of sight. Recently, Frail et 
al. (2001) derived opening angles of GRBs by analysing gra- 
dient changes in the multicolour light curves of GRB after- 
glows. They concluded that the afterglow emission of GRBs 
is generally beamed, with values for 9q ranging between 0.05 
and 0.4 radians. This could hence be an important effect, 
significantly reducing the intensity of the infrared radiation 
emitted from a reprocessed afterglow as compared with the 
predictions in the isotropic model presented here. 

In a preprint pointed out by the referee, Fruchter, Krolik 
& Rhoads (2001) discuss dust destruction by X-rays. They 
found that heating by the optical/UV flash and by X-rays 
have roughly comparable importance for grain evaporation, 
and either may be more important, depending on the burst 
spectrum and grain size distribution. Optical/UV heating is 
generally more important for small grains and X-ray heating 
for large grains, although the division depends on grain com- 
position. They also found that under most circumstances 
grain charging and electrostatic shattering is a more effec- 
tive mechanism in destroying dust grains, compared to heat- 
ing and evaporation, in disagreement with the conclusions 
of Waxman & Draine (2000), who found that grain shatter- 
ing does not appreciably reduce the UV/optical extinction. 
The processes described by Fruchter et al. (2001) could re- 
duce the number of dust grains surrounding a GRB, and so 
further reduce the infrared dust emission. 

If the distribution of the dust is not homogeneous but 
clumpy, then the heating of a dust grain will not only depend 
on R but also on 9 and (f> (F — F[R, a, t, 9, </)]) and equation 
^| will no longer be valid. At some places, radiation levels 
inside the clumpy medium are greatly enhanced with respect 
to a homogeneous medium due to relatively dust-free lines of 
sight, while at other places radiation levels are lower due to 
more effective screening (for example, see Hobson & Padman 
1993; Hobson & Scheuer 1993). 



3.3 Optical depth 

Assume a source of radiation is surrounded by spherical 
shells of dust. The optical/UV flux from the source inci- 
dent on a dust grain at a given radius from the source will 
be attenuated by the dust at smaller radii, and so the opti- 
cal depth, r, between the source and the grain will be non 
zero. In general, t u is defined as 



dr = 



n(a, r)a v (a) da dr 



(10) 



where a v (a) is the absorption cross-section of a dust grain 
with size a and n(a, r) is the number density of dust grains 
with sizes between a and a + da at a distance r from the 
GRB. Because some of the dust grains can be sublimated 
due to intense heating from the source, t u will in general 
also be a function of time. The absorption cross section is 



<7„(a, r, t) 



Q v ,zb B {a)ira 
0, 



t < t dcs (a,r) 
t > *des(o,r), 



(11) 



where tdcs{a,r) is the destruction time of a grain with size 
a at a distance r from the source and <2„, a bs(ffl) the absorp- 
tion coefficient. To compute k v an expression for n(a, r) is 
required. In this model the empirical law 



n(a, r) = C(r)a 



(12) 



(Mathis, Rumpl & Nordsieck 1977; Hildebrand 1983) was 
used for the grain size distribution, in which 



C(r) = 



3 n H (r) m H 



8npM g /M d f^dmax - \/a min ) 



(13) 



nn(r) is the number density of neutral hydrogen, mu is the 
atomic mass of hydrogen, M s /M d is the gas-to-dust mass 
ratio, p is the mean mass density of the dust grains, and 
cimin and a max are the minimum and maximum grain radii. 

Using equations (|l^) to (^), the total optical depth 
can now be calculated. This is the optical depth due to dust 
immediately surrounding the GRB. In addition, an after- 
glow light curve can be attenuated by dust along the line of 
sight within the host galaxy of the burst, in the intergalactic 
medium (Perna & Aguirre 2000) and in our own galaxy. 

With the expressions for the light curve of GRB after- 
glows with or without an optical flash (see section | above) , 
the luminosity in the photon energy range hv = 1 — 7.5 eV, 
Li_7.5(t), can now be calculated, and so the dust proper- 
ties T(t), r(r,t) and a(t) can be computed throughout the 
enshrouding dust cloud, as the optical emission of the GRB 
travels through it. 



3.4 Numerical model 

3.4-1 Grain material 

For the simulations a model of dust consisting of carbon 
grains with various sizes was used. Graphite grains are bet- 
ter absorbers of optical and UV photons than silicate grains 
(Mezger, Mathis & Panagia 1982; Draine & Lee 1984). On 
the other hand, silicate grains are destroyed at lower tem- 
peratures than graphite grains because their chemical bonds 
are weaker (Guhathakurta & Draine 1989; Evans 1993). 
Graphite grains will thus get hotter, but are also more ro- 
bust, and so the destruction time of graphites and silicates 
are not expected to differ by large factors. Hence, because of 
their higher absorptivity, graphite grains are likely to dom- 
inate the optical depth at optical and UV wavelengths. 
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3.4.2 Absorption and emission coefficients 

The value of the absorption coefficient of a graphite grain 
in the energy range f — 7.5 eV, Q a b s , depends on the grain 
size. Based on the results of computations by Draine & Lee 
(f984), if we define an effective wavelength A c ff to account 
for the interaction of the physical properties of the grains 
and the illuminating radiation field, then the absorption co- 
efficient 



Qabs (a) 



1 U a A efr , a 



< O.lAefi 
a > 0.1 A c ff. 



(14) 



A value of 0.5 /im is taken for A c ft , because in the GRB 
afterglow model approximately half the energy is radiated 
at longer wavelengths and half at shorter wavelengths, de- 
pending on whether or not there is an optical flash, which 
is expected to be bluer than a typical afterglow. Changing 
A e ff from 0.4 /im to 0.6 /im made little difference to the out- 
come of the simulations: the destruction time of the grains 
and the dust emission both changed by only about f per 
cent. Empirical values from Draine (200f) were taken for 
the emission coefficients, Q cm (equation^). 

3-4-3 Density profiles and dust mass 

Two different dust density profiles were considered; a con- 
stant density at all radii and an isothermal density profile, 
that is n(r) oc r~ 2 . A gas-to-dust mass ratio of fO 2 was as- 
sumed (Hildebrand f983; Evans f993). In the model it is 
assumed there is no dust within 0.3 pc of the GRB. This is 
the Stromgren radius of an HII region with riH = 5000 cm -3 
surrounding a star of spectral type 06. If the progenitors 
of GRBs are massive stars of this type, then it is unlikely 
that significant amounts of dust will be able to survive 
within the Stromgren sphere. Assuming a different inner 
radius, however, does not affect the optical depth or the 
dust mass strongly. Even if the inner radius is reduced to 
zero, then the dust mass will change by only a few percent. 
The maximum radius is taken to be 3.67 pc and 10.33 pc in 
the case of a constant and isothermal density profile respec- 
tively. The total dust mass is 0.234 (nn 3 /cm -3 ) Mq and 
0.016 (nH 3 /cm -3 ) Mq in these cases respectively, where 
the density of hydrogen at 0.3 pc (wh 3 ) is a free parameter 
in the model. 



3.4-4 Minimum and maximum grain size 
As mentioned in section 



3.3 an empirical law (equation 



|l2] ) was used for the distribution of dust grain sizes from 
0.005 /im (a m i n ) to 2.5 fim (a max ). The minimum grain size 
was derived by Mathis et al. (1977). The maximum grain size 
is adopted from Witt, Oliveri & Schild (1990) and Aitken 
et al. (1993). However, the maximum grain size assumed by 
different authors varies widely, in papers about theory as 
well as in papers about observations. To investigate the ef- 
fect of the minimum and maximum value of the grain size 
on the results of the dust model, the range was extended to 
0.001 /im as a lower limit, to reach molecular sizes, and to 
5/im as an upper limit. 

The upper limit to the grain size has the strongest influ- 
ence on the total number of dust grains for a fixed dust mass 



Table 1. Optical depth all the way through the dust cloud 
for different profiles and limits of grain sizes. The total dust 
mass and maximum radius in the case of a flat density profile 
is 0.234 (n,[j 3 /cm -3 ) Mq and 3.67pc. In the case of decreas- 
ing profile these are 0.016 (njj 3 /cm -3 ) Mq and 10.33 pc. To 
calculate the optical depth at other wavelengths, it is assumed 
that the optical depth varies as oc A -1 at wavelengths 0.5 (im 
< A < 1.25 (im (for example, Evans 1993) and as oc A -16 at wave- 
lengths A > 1.25 (im (Cardilli, Clayton &; Mathis 1989; Mathis 
1990). 
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in the model, because the maximum grain size determines 
the value of C(r) (equation [l3|). If the maximum size of the 
grains is increased by a factor of 2, then the value of C(r) will 
decrease by a factor of approximately y2. This influences 
both the optical depth and the dust luminosity; for exam- 
ple, for the same total dust mass, the decrease in the number 
of small grains as a max increases will lead to a reduction in 
the optical depth, because the optical depth is dominated 
by small er g rains that are still larger than molecules (see 
equation pdf dr ~ a -1 ' 5 da), as shown in Table|l[ Hence, 
the choice of maximum grain size has a significant effect, 
because it determines the number of dust grains. However, 
if C(r), and thus the number of grains with a given radius, is 
kept fixed, then neither the optical depth nor the dust emis- 
sion change significantly as a function of a max : the change in 
dust emission with a max is less than 1 per cent at all wave- 
lengths. The effect of changing the upper limit of grain size 
in the model is thus negligible if the normalization of n(a, r) 
is kept constant, although the total mass of dust will then 
change. 

Decreasing the lower limit of the grain size from 
0.005 /im to 0.001 /im with 71h 3 fixed increases the opti- 
cal depth by about 11 per cent (see Table Q), while C(r) 
only changes by about 2 per cent. Because of the increased 
amount of energy absorbed by the smallest grains when a m i n 
is reduced, the optical/UV flux density of the late-time af- 
terglow is also reduced by 40 per cent. However, the dust 
luminosity does not change by a large amount, because the 
fractional change in the amount of energy absorbed is only 
of order 0.1 per cent. Decreasing the minimum grain size fur- 
ther has a negligible effect on both the dust emission and the 
optical depth, as the additional amount of energy absorbed 
is very small. 

Absorption by large molecules or Polycyclic Aromatic 
Hydrocarbon maolecules (PAHs) only become an important 
factor for the optical depth at photon energies hv > 10 eV, 
at which their absorption cross sections increase to be com- 
parable to their physical cross section (Draine & Lee 1984). 
Below those photon energies, the optical depth is dominated 
by bigger grains (see also Desert, Boulanger & Puget 1990). 
Because photons with energy hv > 7.5 eV will mainly be 
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Table 2. Values of some important parameters, used for simu- 
lations. In the table OF is the representation of the optical flash 
(e.g. Gaussian), L pea k is the peak luminosity of the optical flash 
in the energy band 1 — 7.5eV in ergs -1 and profile is the density 
profile. The duration of the optical flash is in each case 10 s. In 
Fig. [l] the unattenuated light curves are plotted. 



light curve 



Gaussian OF 
afterglow 



Gaussian OF 
afterglow 



afterglow only 



Top hat OF 
afterglow 





riu /cm 3 


profile 


name 




5xl0 3 


flat 


Ai 




5xl0 4 


flat 
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5xl0 3 


isothermal 


A3 




5xl0 4 


isothermal 


A 4 




5xl0 5 


isothermal 


A 5 




5xl0 3 


flat 


Bi 




5xl0 4 


flat 


B 2 


10 48 


5xl0 3 


isothermal 


B 3 




5xl0 4 


isothermal 


Ba 




5xl0 5 


isothermal 


Br, 




5xl0 3 


flat 


Ci 




5xl0 4 


flat 


C-2 




5xl0 3 


isothermal 


c 3 




5xl0 4 


isothermal 


r* 4 




5xl0 5 


isothermal 


C 5 


10 48 


5xl0 5 


isothermal 


D 



absorbed both by H and H2 in the early phases of the GRB 
(Waxman & Draine 2000, see section [|), the smallest grains 
and molecules are unlikely to contribute significantly to the 
total dust emission. 

The minimum and maximum grain sizes used in the 
model were assumed to be 0.005 /im and 2.5 /im respectively. 
However, changing the minimum and maximum grain size 
was found not to influence the results of the model signifi- 
cantly, if the number of grains with a given radius was kept 
fixed. 



4 RESULTS: OPTICAL LIGHT CURVES 

In Table ^ the parameters describing the light curves and the 
properties of surrounding dust clouds used here are listed. In 
Fig.|i|a light curves are plotted for a GRB with an afterglow 
and an optical flash with a peak luminosity of 10 49 ergs _1 
as it would appear to an observer at a distance of 1 Mpc, 
including the effects of absorption by dust surrounding the 
GRB. 

As can be seen in Fig. the optical flash is able to de- 
stroy large amounts of dust. Although the total dust mass 
in the clou d varies a great deal from model to model (see 
section 3.4), and can be as great as 2300 Mq, the visibility 
of the afterglow is hardly reduced after the first few seconds. 
In the case of a flat dust density profile with tih = 0.5 and 



5x10 cm" , as represented by curve A\ and A2 respectively, 
all the dust is destroyed by the optical flash and observations 
of the afterglow would not reveal any dust surrounding the 
GRB. In the case of an isothermal density profile, some of 
the bigger dust grains in the outer shells survive the passage 
of the flash. This is not because of screening by the absorp- 
tion due to grains at smaller radii, but because the distance 
to the burst is sufficiently large that the intensity of the op- 
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Figure 2. (a) Integrated flux of a GRB afterglow plus a Gaussian 
optical flash with L peak = 10 49 erg s 



between v = 2.4 x 10 14 Hz and v 



1 at a distance of 1 Mpc 
..8 X 10 15 Hz. The dotted 
line is the light curve with no extinction. The thin solid line and 
the thick solid line include absorption by dust with a flat density 
profile and a density of njj = 0.5 and 5x 10 4 cm~ 3 respectively 
(light curve A\ and A2 of Table|i|). The dash-dotted line, the 
dashed line and the long-dashed line include absorption by dust 
with an isothermal density profile and a maximum density of 



"H,n 



0.5, 5 and 50xl0 4 cm 3 respectively (light curve A3, 



Ai and A5 of Tablebl). (b) As (a), but with an optical flash with 
Lpeak = 10 48 ergs — 1 . The light curve which includes absorption 
by dust with constant density of njj = 5 X 10 4 cm~ 3 (light curve 
B2) is not shown, because the flux is close to zero because of 
the very powerful absorption in this case, (c) As (a), but with no 
optical flash. Light curves C2 and C5 are not shown, because the 
flux is almost zero due to absorption by the dust. 
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tical/UV flux has decreased below the level at which grains 
are destroyed by sublimation. For example, beyond about 
4.5 pc the flux density of even an unattenuated flash with a 
peak luminosity of 10 49 ergs _1 is not sufficient to sublime 
the largest grains (a > 1 /im). 

When the peak of the optical flash is less intense, 
10 48 ergs _1 , less dust is destroyed along the line of sight, 
and so the afterglow is attenuated to a greater degree (see 
Fig.|^b). In the case of a flat profile with tih = 5x 1 4 cm -3 , 
practically all the optical/UV energy of the flash and after- 
glow is absorbed by the surrounding dust. In the case of an 
isothermal profile and a density of wh = 5x 10 5 cm~ 3 , the 
afterglow is attenuated by a factor of about 2.5 x I0 5 or If 
magnitudes. Hence, it is very unlikely that an optical/UV 
afterglow would be detected; for example, the unattenuated 
R-band magnitude of a GRB afterglow at a relatively low 
redshift z = 0.3 one day after the start of the burst, which 
is predicted to be R ~ 18, would be reduced to R ~ 29. 

In the densest environments, the afterglow of a GRB 
with a relatively weak optical flash will be heavily attenu- 
ated; however, in less dense environments, such a burst will 
only be marginally affected. 

If there is no optical flash (or if the intensity of the op- 
tical flash has an insignificant value of I/ pea k ;$ 10 47 ergs -1 ), 
then almost no dust is destroyed, and so the optical/UV af- 
terglow emission from the GRB will be severely attenuated 
in all cases: see Fig. |2|c. If the shape of the optical flash is 
a top-hat and not a Gaussian, then more dust will be de- 
stroyed, simply because there is more energy in a flash with 
a top-hat profile than in a Gaussian flash with the same 
peak flux. 

In conclusion, the higher the peak luminosity of the op- 
tical flash, the greater is the amount of dust that is expected 
to be destroyed in the environment of the GRB, and so the 
afterglow will be less attenuated. If the peak luminosity of 
the optical flash is > I0 49 ergs -1 , then the intensity of the 
afterglow emission will never be reduced significantly. If the 
peak luminosity of the optical flash < 10 47 ergs _1 , then the 
afterglow will always be highly obscured, unless the total 
mass of the dust within a few pc surrounding the GRB is less 
than IOOM0, which would be the case for tih < 1 3 cm -3 
and a gas-to-dust mass ratio of about I0 2 . 



5 RESULTS: DUST EMISSION 

5.1 Dust emission from the GRB environment 

The dust emission observed from around a GRB at a given 
wavelength depends on the density and profile of the sur- 
rounding dust, the peak luminosity of the optical flash, the 
light curve of the GRB afterglow and the redshift of the 
GRB. To illustrate the key properties of the dust emission, 
a single optimistic set of model parameters was assumed. 
The density profile of the dust was assumed to be flat with 
nu = 10 4 cm -3 . The GRB optical/UV light curve was as- 
sumed to be i?2 in Table ^| The redshift of the GRB was 
assumed to be z = I, corresponding to a luminosity distance 
d L = 5.4Gpc (H = eSkms^Mpc" 1 and g = 0.5). The 
GRB emission was assumed to be isotropic, not beamed. 
This optimistic set of parameters also ensures an intense 
dust emission signal - more than 2400 Mq of dust surrounds 
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Figure 3. (a) Emission from dust surrounding a GRB at z = 1. 
Light curve E>2, which is described in Table 0, was assumed to 
heat the dust, which absorbed all the energy of the light curve 
as shown in Fig. Here the reprocessed emission at observed 
wavelengths of 1.6, 2.2, 3.5, 4.8 and 10 fim are plotted, (b) The 
corresponding emission at four longer observed wavelengths of 20, 
60, 90 and 170 um. In both figures the flux densities at the same 
wavelengths of the luminous dusty star-forming galaxy Arp 220, 
redshifted to z = 1, are shown for comparison. 



the GRB, and nearly all the energy of the afterglow and opti- 
cal flash is absorbed and re-emitted by the dust. In Fig. ^ the 
associated dust emission predicted by the radiative transfer 
model at eight observed wavelengths between 1.6 /zm and 
170 jUm is plotted. A characteristic feature of the emission 
from the spherical distribution of dust in this case is that 
the detected signal is a top-hat function. The emission at the 
high-energy end of the dust radiation spectrum is radiated 
predominantly by hot grains near their sublimation tem- 
perature in the time before they vaporize. These grains are 
located close to the GRB, at a distance R ~ 1 pc. Because 
these grains are vaporized by the intense radiation of the op- 
tical flash, the dust emission at short wavelengths starts to 
decrease about t > 2 (1 + z) R/c « 6(1 + z) years after the 
GRB is detected. The emission at longer wavelengths is pre- 
dominantly due to cooler dust, which is further away from 
the GRB. The dust cloud in this model was assumed to have 
a maximum radius i? ma x = 3.67 pc, and so the emission de- 
creases rapidly after t > 2 (1 + z) i? m ax/c ~ 22 (1 + z) years. 
The dust emission spectrum is predicted to peak at a rest- 
frame wavelength of about 8 /jm. 
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Figure 4. Emission from dust surrounding a GRB at z = 1, in- 
cluding emission from the GRB afterglow at five observed wave- 
lengths. The light curve that heated the dust is light curve B2 as 
shown in Fig. I (see Table|| for its characteristics). 

The dust emission predicted at the observed wavelength 
170 fim at t < 50 days is due to the dust surrounding the 
GRB radiating at its original temperature in the interstellar 
radiation field of the GRB host galaxy, which is assumed to 
be 30 K. If the dust temperature assumed is different, then 
the flux changes by only a small factor at this wavelength. 
At rest-frame wavelengths longer than 100 /im, at which the 
spectral energy distribution of dust at 30 K peaks, there is 
very little excess emission due to additional heating by the 
GRB. 

For comparison, the expected flux densities of an 
ULIRG similar to Arp 220 (Klaas et al. 1997), at a redshift 
z = 1, is also plotted in Fig. ^j. The intensity of the emis- 
sion from dust heated by the GRB at rest-frame wavelengths 
shorter than about 2.5 /xm is expected to be less than that 
from starlight in such a galaxy. However, the dust emission 
from the GRB at wavelengths of about 5 fim can be an order 
of magnitude greater than the emission from the galaxy. At 
rest-frame wavelengths longer than about 25 ^m, infrared 
emission from dust heated by star formation at tempera- 
tures T ~ 40 K rises to dominate the spectrum. Therefore, 
the dust emission from a GRB within a galaxy like Arp 
220 at 2 = 1 could best be observed at mid-infrared wave- 
lengths between 5 /jm and 20 ^m. Note, however, that Arp 
220 is a very luminous galaxy. It is likely that the host galaxy 
of a GRB will typically be less luminous, which would im- 
prove the observability of the GRB dust emission. Also, the 
infrared emission of a galaxy will be much more extended 
than emission from the surroundings of a GRB and could 
be distinguished if high spatial resolution observations were 
available, such as would be possible using either the Ata- 
cama Large Millemeter Array (ALMA) or a space-borne far- 
infrared interferometer like SPECS (see section ^ for more 
details). 

5.2 The infrared GRB synchrotron afterglow 

In Fig.^J only dust emission from the environment of the 
GRB was plotted. However, it is likely that in the first few 
days after the GRB, the observed spectrum at mid-infrared 



wavelengths will be dominated by the GRB afterglow. In 
Fig. ^ the dust emission at five different wavelengths is plot- 
ted, including the emission from the GRB afterglow. The 
light curves of the GRB dust emission are the same as plot- 
ted in Fig.|[ It is assumed that the infrared afterglow is not 
attenuated by dust surrounding the GRB, however, at the 
shortest infrared wavelengths, dust along the line of sight 
could reduce the flux density of the afterglow significantly, 
as shown in Table jjj. After a few days, the afterglow emis- 
sion is decreasing as t~ 1/i75 (Sari et al. 1998). All the light 
curves shown in Fig. ^ are dominated by the afterglow in the 
first few weeks after the start of the GRB. After the first 
few weeks, the light curve flattens as the afterglow emis- 
sion decreases below the level of the heated dust emission, 
which then dominates the light curve at late times. If the 
flux density of the dust emission is ten times lower than that 
assumed in this model, then the afterglow will dominate the 
spectrum for if) 1 / 1,375 « 5 times longer. In other less opti- 
mistic models the GRB dust emission is typically lower and 
so the afterglow will dominate the spectrum for a longer 
period. 

5.3 Dust emission in various environments 

The results above showed that the visibility of the optical 
afterglow can depend strongly on the GRB environment. 
Different amounts of energy are absorbed by different dis- 
tributions of dust, thus influencing the infrared emission. In 
Fig.|B|a the peak flux of the dust emission spectrum is calcu- 
lated for the ^4-series models listed in Table ^ at z = 1. The 
maximum flux is expected to occur roughly (1+z) years after 
the GRB is detected. The time profile of the dust emission is 
similar to that displayed in Fi g.pt The spectrum of the dust 
emission described in section |5.l| is plotted for comparison. 
In the case of flat dust density profiles, almost all the dust 
is vaporized, and so the spectrum is dominated by emission 
from dust near the sublimation temperature, peaking at a 
wavelength about (1 + z) /im. In the other cases, the ob- 
served emission peaks at a longer observed wavelength close 
to 10 (1 + z) fim. 

Because the optical flash is expected to destroy a large 
amount of dust in the A-series light curve models, the in- 
tensity of the dust emission is less than that predicted in a 
model with a less luminous optical flash. This is illustrated 
in Fig. ^3, in which the dust emission spectrum is plotted 
for dust heated by the same GRB afterglow, but with an 
optical flash that is ten times weaker, represented by the 
B-series models, see Table ^. Light curve D has the same 
parameters as light curve B5, except that the optical flash 
is this time described by a top-hat. Since in the model the 
total energy radiated in the optical flash is greater, more 
dust is destroyed and so the flux density of the dust emis- 
sion is smaller. Because less dust is destroyed in the B-series 
models as compared with the A-series models, the intensity 
of the dust emission is greater in the B-series models for a 
given density profile. 

The emission is expected to be even greater if the dust 
is heated by a GRB with an afterglow but no optical flash, 
as shown for the C-series models in Fig.^c. If the peak lu- 
minosity of the optical flash is small, then almost no dust 
is destroyed and so the dust emission is intense. Except for 
models Ai and A2, for which all the dust is destroyed by the 
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Figure 5. (a) The peak intensity of the emission spectrum of dust heated by an optical flash with Z/ pca k = 10 49 ergs -1 and an afterglow 
of a GRB at Z = 1. The dust parameters are described in Table For comparison the spectrum of the more powerful dust emission in 
model B2 is also plotted, (b) The peak intensity of the emission spectrum of dust heated by an optical flash with Z/ pe ak 

= 10 48 erg s 1 

and an afterglow of a GRB at z = 1 in different dusty environments, described in Table ^. Curve D has the same parameters as curve 
B5, except that the optical flash is not a Gaussian but a top-hat. (c) The peak intensity of the emission spectrum of dust heated by an 
afterglow of a GRB at z = 1 with no optical flash, described by parameters C\ to C5 in Table g. For comparison the spectrum of the 
dust emission in model B2 is also plotted. 



optical flash, dust with a flat density profile emits more ra- 
diation per unit mass than dust with an isothermal profile. 
To obtain the greatest dust emission from a GRB, an optical 
flash with a low peak luminosity, a high-density environment 
and a flat density profile are all required. 



5.4 Dust emission and the redshift of the GRB 

In Fig. ^| the predicted maximum intensity of emission by 
dust heated by a GRB, with a model B2 optical flash after- 
glow, are plotted in various bands as a function of redshift. 
In this optimistic model there is strong evolution with red- 
shift: the predicted flux density decreases rapidly as the red- 
shift increases, especially at the short wavelengths. At longer 
wavelengths (A > 20 (im) the flux density decreases more 
slowly, as the observed wavelength approaches the peak of 
rest-frame spectrum. 



6 OBSERVABILITY OF DUST EMISSION 
FROM GRBS 

The observed wavelengths 3.5, 4.8, 5.8 and 8.0 /im in Fig.^, 
correspond to bandpasses of the InfraRed Array Cam- 
era (IRAC) on board the Space Infrared Telescope Facility 
(SIRTF^\). The Multiband Imaging Photometer for SIRTF 
(MIPS) has bandpasses at 24, 70 and 160 fim. In 1 hour 



of integration time, SIRTF could detect the near-infrared 
dust emission from a GRB at a redshift 2 < 2. At the 
longer wavelengths observable with the MIPS instrument, 
it is unlikely that SIRTF will detect dust emission from 
a GRB at a redshift z > 1, because the emission is too 
weak. Also, at the longest wavelengths observable for the 
MIPS instrument, 70 and 160 /im, confusion noise from un- 
resolved distant galaxies will increase the minimum detec- 
tion depth, reducing the maximum redshift at which dust 
emission could be observable (see Blain 1999). Because the 
GRB dust emission at far-infrared and submillimetre wave- 
lengths is expected to be much less intense than that at 
near-infrared wavelengths, facilities operating in those wave- 
bands, like the future Stratospheric Observatory for Infrared 
Astronomy (SOFIAjjj) and the Far InfraRed and Submillime- 
tre Telescope (FIRST-Herschel^, are likely to be unsuc- 
cessful in detecting dust emission from high redshift GRBs. 
However, it should be possible to resolve the GRB dust emis- 
sion using high spatial resolution facilities like ALMA^] and 
the SPECS (see Leisawitz et al. 1999, 2000). In cases where 
the optical afterglow is heavily obscured by dust and thus in- 
visible to optical searches, ALMA could provide a very accu- 
rate position for the GRB afterglow, once the GRB position 
is known to within a few arcseconds. Accurate positions of 
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The homepage 



Up-to-date in form at 



http:/ /sirtf.caltech.edu/ 



about SIRTF can be found at 



general/Proj ects/First / first .html 
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Figure 6. (a) The maximum flux density of dust emission from a 
GRB as a function of redshift at several different observed wave- 
lengths model £?2 ■ The flux density decreases rapidly as the red- 
shift increases, between redshifts 1 and 2 the flux density de- 
creases by a factor of about ten. The top four curves are plotted 
for wavelengths observable using the SIRTF-TRAC camera. The 
maximum redshifts at which the IRAC camera could detect the 
dust emission with a signal to noise ratio of 5 within 1 hour of 
integration time are 1.7, 2.1, 1.8 and 2.0 for the 3.6, 4.5, 5.8 and 
8.0 bandpass respectively, (b) The maximum flux density of dust 
emission from the same GRB at several longer observed wave- 
lengths. The MIPS instrument on SIRTF has three observing 
bands, at 24, 70 and 160 ^m. The maximum redshifts at which 
MIPS could detect the dust emission with a signal-to-noise ra- 
tio of 5 within 1 hour of integration time are 0.8, 0.4 and 0.1 for 
the 24, 70 and 160 /im bandpass respectively Confusion noise will 
limit the depth of observations at the longest wavelengths, 70 and 
160 /^m (Blain 1999 and references therein). 



optically obscured afterglows are currently provided by the 
VLA at radio frequencies, as shown by Frail et al. (2000). 
Such accurate positions can be provided by the Swift satel- 
which will locate a GRB precisely at X-ray energies 
within 2.5 arcseconds (Gehrels et al. 1999). Already, a fad- 
ing counterpart to a GRB with a flux density of several mjy 
has been detected at 850 pm (Smith et al. 1999). In a ten- 
minute observation ALMA will have a 5cr detection limit of 
about 0.2 mjy at a wavelength of 870 /im. This will be deep 



II The homenag-e of Swift w ith the latest updates can be found 
at http://swift.sonoma.edu/ 



enough to detect the submillimetre-wave afterglow of a GRB 
as far away as a redshift of 10 within a day of the start of 
the GRB, while the afterglow of a GRB at a redshift of 5 
would be detectable for three weeks. 



7 DISCUSSION 

In this paper the visibility of optical/UV GRB emission ab- 
sorbed by dust and the intensities of the reprocessed radi- 
ation emitted by dust were predicted. The heating effect of 
GRB afterglow and prompt optical flash was calculated in 
a model whose parameters are consistent with the relatively 
sparse existing observations. A radiative transfer model was 
used to simulate the interaction between GRB radiation and 
dust grains; the temperature, opacity and sublimation rate 
of the dust grains in the environment of a GRB, and the 
infrared emission from this dust was calculated. Note that 
in the model the GRB emission was assumed to be isotropic 
and that dust destruction by X-rays was neglected. If GRBs 
are beamed and/or grains are predominantly destroyed by 
X-rays, then the dust emission as computed in the model 
could be reduced significantly. 

The simulations showed that the visibility of the opti- 
cal flash and afterglow of a GRB depends strongly on the 
characteristics of both the surrounding dust and the GRB 
emission itself. If a GRB is surrounded by a large amount 
of dust (Afdust > 100 Mq), then it is unlikely that an opti- 
cal afterglow will be observed, unless the GRB is accompa- 
nied by a very powerful optical flash with a peak luminos- 
ity Lpeak > 10 49 ergs^ 1 , which would destroy almost all the 
dust along the line of sight within a few seconds of the GRB. 
Because of the large-scale destruction of dust grains in this 
scenario, the intensity of the associated dust emission would 
be quite weak. If the optical flash associated with a GRB is 
less powerful, then less dust is destroyed. The opacity is thus 
higher, and so it is possible that the energy in the afterglow 
is almost completely absorbed, and that no late-time op- 
tical emission from the GRB would be observed. However, 
the reradiated energy might be detectable. In particular, 
the IRAC camera on SIRTF could detect the near-infrared 
dust emission in such cases out to redshifts of about 2, if 
the GRB could be located to within its 5— arcmin field of 
view. Because the light curve of this dust emission does not 
have any specific temporal structures, and could therefore 
be confused with the flux density of dust radiation from the 
general interstellar medium of the host galaxy, an accurate 
position and flux density of the GRB afterglow would need 
to be provided by other observations. At present these can 
be carried out with the VLA. In the future high-resolution 
submillimetre-wave observations with ALMA will be able to 
confirm that the origin of the infrared energy was a GRB. 

A combination of observations is needed to test whether 
a GRB is surrounded by a large amount of dust. First of all, 
the lack of a detection of an afterglow in the optical wave- 
band could be a first indication that there is a large mass of 
dust near the GRB. If radio observations provide an accu- 
rate position for the afterglow, then the SIRTF mid-infrared 
telescope could be used to try to detect the reprocessed dust 
emission directly. However, only in the most optimistic sce- 
narios, in which more than 10 3 Mq of dust is present and the 
optical flash of the GRB is weak, is infrared emission likely 
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to be observed from dust heated by the optical flash and af- 
terglow of the GRB. If either the GRB emission is beamed, 
or X-ray grain destruction is important, or less than 10 2 Mq 
of dust surrounds the GRB, then the infrared dust emission 
is likely to fall below the detection limits of SIRTF unless 
the GRB is at a low redshift z < 0.5. 

The results of the model can be compared with the 
results of existing observations for optical flashes from 
GRBs. A powerful optical flash was observed to accom- 
pany GRB990123, detected twenty seconds after the start 
of the GRB (Akerlof et al. 1999). From the observed red- 
shift z = l.G (Galama et al. 1999b), the peak luminosity 
was calculated to be ~ 10 50 ergs _1 . This flash could easily 
have vaporized more than 10 3 Mq of dust around the GRB 
within a few seconds. The late-time afterglow for this GRB 
is thus unlikely to have been attenuated by dust surrounding 
the GRB, and so the dust emission in this case was prob- 
ably very low. Attempts to detect the prompt optical flash 
from other GRBs have provided only upper limits (Akerlof 
et al. 2000). This does not necessarily imply that the opti- 
cal emission in those cases was extinguished by surrounding 
dust. However, if the peak luminosity of the optical flash of 
a GRB is well correlated with either gamma-ray fluence or 
gamma-ray peak flux, then in two of the observed cases the 
optical flash would have I/ pea k — 10 49 ergs -1 . This luminos- 
ity is great enough to sublime most of the surrounding dust, 
and thus lead to the detection of optical emission. Even if 
these GRBs occurred in very dusty environments, then the 
lack of optical emission indicates that the peak luminosity 
of the optical flash is probably not strongly correlated with 
either gamma-ray fluence or gamma-ray peak flux. 
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